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Abstract 

We present a phenomenological study of r-sleptons fi,2 and r-sneutrinos i/r in the Minimal Su- 
persymmetric Standard Model with complex parameters At, jJ. and Mi. We analyse production 
and decays of the fi,2 and Dt at a future e^e~ collider. We present numerical predictions for the 
important decay rates, paying particular attention to their dependence on the complex parameters. 
The branching ratios of the fermionic decays of fi and show a significant phase dependence for 
tan/3^10. For tan/3^10 the branching ratios for the f2 decays into Higgs bosons depend very sen- 
sitively on the phases. We show how information on the phase ipA^ and the other fundamental fi 
parameters can be obtained from measurements of the fi masses, polarized cross sections and bosonic 
and fermionic decay branching ratios, for small and large tan /3 values. We estimate the expected 
errors of these parameters. Given favorable conditions, the error of At is about 10% to 20%, while 
the errors of the remaining stau parameters are in the range of approximately 1% to 3%. We also 
show that the induced electric dipole moment of the r-lepton is well below the current experimental 
limit. 



1 Introduction 



So far most phenomenological studies on supersymmetric (SUSY) particle searches have been performed 
within the Minimal Supersymmetric Standard Model (MSSM) with real SUSY parameters. In this paper 
wc study the production and decays of r-slcptons and r-sncutrinos at an e~^e~ linear collider in the 
MSSM with complex SUSY parameters. 

In the SUSY extension of the Standard Model (SM) one introduces scalar leptons £l, ^r, scalar 
neutrinos £>£ and scalar quarks ql, Qr. as the SUSY partners of the leptons £l,r, neutrinos vi and quarks 
qL,R, respectively For each definite fcrmion flavor the states /l and f^ are mixed by Yukawa terms. 
The mass eigenstates are /i and /2, with m^^ < mj^ P]. For the sfermions of the first and second 
generation — mixing can be neglected. For the third generation sfermions, however, — //} mixing 
has to be taken into account due to the larger Yukawa coupling |31 ^ . 

In the case of the r-sleptons — fji mixing is important if the SUSY parameter tan/3 is large, 
tan/3 > 20. The lower mass eigenvalue rrif^ can be rather small and the fi could be the lightest charged 
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SUSY particle. The experimental search for the r-sleptons and the r-sncutrino and the determination of 
their parameters is, therefore, an important issue at all present and future colliders. Pair production of 
r-sleptons and r-sneutrinos will be particularly interesting at an e'^e~ linear collider with centre of mass 
energy -y/s = 0.5 — 1.2 TeV. At such a colHder and with an integrated luminosity of about 500 fb^^ it 
will be possible to measure masses, cross sections and decay branching ratios with high precision [SJIS]. 
This will allow us to obtain information on the fundamental soft SUSY breaking parameters of the third 
generation slepton system. 

In the recent phenomenological study of 3rd generation sfermions in the real MSSM it has been 
shown how the masses and the mixing angle of the stop system can be determined by measurements 
of the production cross sections with polarized beams The results of a simulation of e~^e~ — > titi 
with the decay modes ii — > x^c and ii — > xtb and including full SM background in jSj imply that with 
an integrated luminosity of 500 fb^^ an accuracy of the order of 1% or better may be obtained. The 
numerical precision to be expected for the determination of the underlying SUSY parameters Mq, Mjy 
and (real) At has also been given. For low tan /3 one can expect similar results for the sbottom and stau 
systems [SlISIZl- 

The assumption of real SUSY parameters has partly been justified by the very small experimental 
upper limits on the electric dipole moments (EDM) of electron and neutron. A possibility to avoid the 
EDM constraints is to assume that the masses of the first and second generation sfermions are large 
(above the TeV scale), while the masses of the third generation sfermions are small (below 1 TeV) 
Pj. Another possibility is suggested by recent analyses of the EDMs, which have shown that strong 
cancellations between the different SUSY contributions to the EDMs can occur ^Hl ■ As a consequence of 
these cancellations it has turned out that the complex phase of the Higgs-higgsino mass parameter jj, is 
much less restricted than previously assumed, whereas the complex phases of the soft-breaking trilinear 
scalar coupling parameters Af are practically unconstrained I12j . For example, in a mSUGRA- 
type model with universal parameters M1/2, Mq, tan/3 and complex Aq, with being determined by 
radiative electroweak symmetry breaking, the phase of /i is constrained to |(/3^|<0.1 — 0.2 for low values 
of the scalar mass parameter, Mo<AOO GeV, and becomes less constrained for higher values of Mq. The 
phase of Aq, (pA(n turns out to be correlated with ip^, but otherwise not restricted |12l In models 
with more general parameter specifications also ip^ turns out to be less constrained |14| . In any case, 
this means that in a complete phenomenological analysis of production and decays of third generation 
sfermions one has to take into account that the SUSY parameters fi and Af may be complex and one 
has to study the implications that follow for the important observables. 

In our present phenomenological study of 3rd generation sleptons we use the MSSM as general frame- 
work and we assume that the parameters /x, At and Mi are complex (A-r is the trilinear scalar coupling 
parameter of the f^-system and Mi is the U{1) gaugino mass parameter). We neglect flavor changing 
CP violating phases and assume that the scalar mass matrices and trilinear scalar coupling parameters 
are flavor diagonal. We perform an analysis of production and decay rates of fi, f2 and i^t at an e+e~ 
linear collider with a CMS energy y/s = 0.5 — 1.2 TeV. We include also explicit CP violation in the Higgs 
sector induced by stop and sbottom loops with complex parameters as in E| and jE] , using the 
loop-corrected formulae of ^Hl- Our present study is an extension of the corresponding one in the MSSM 
with real parameters in [3. Compared to the real MSSM, the inclusion of the complex phases (fA^, 
and fu{i) of ^t, ^J■ and Mi means that the number of independent fundamental SUSY parameters is 
increased. In order to determine all these parameters one has to measure more independent observables 
than in the real case. 

In principle, the imaginary parts of the complex parameters involved could most directly and unam- 
biguously be determined by measuring suitable CP violating observables. However, in the fi-system this 
is not straightforward, because the are spinless and their main decay modes are two-body decays. A 
possible method has been proposed in |18| . which is applicable if the mass splitting between the mass 
eigenstates f 1 and f2 is very small. If ^ — nif^ is of the order of the decay widths, fi — f2 oscillations will 
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occur which can lead to large CP violating asymmetries in e"'"e~ annihilation. In Ref. an analysis of 
^JL'^ fifj with longitudinally and transversely polarized beams has been given and the observables 
sensitive to CP violation in the Ti sector and Higgs sector have been classified. 

On the other hand, also the CP conserving observables depend on the phases of the underlying complex 
parameters, because the mass eigenvalues and the couplings involved are functions of these parameters. 
In particular, the various decay branching ratios depend in a characteristic way on the complex phases. 
The main purpose of the present paper is a detailed study of the fermionic decay branching ratios of 
fi, f2 and Vr, and the bosonic decay branching ratios of T2 and and their dependences on the phases 
Va^ , and (pu(i) ■ In |20| we have published first results of our study. In the present paper we give the 
analytic expressions for the various decay widths with complex couplings. We present a more detailed 
numerical study of the phase dependences of the various branching ratios. We also discuss how these 
phase dependences can be qualitatively understood on the basis of the analytic expressions for the decay 
widths. Furthermore, we give a theoretical estimate of the precision to be expected for the determination 
of the complex phases together with the other fundamental parameters of the fi-system by measurements 
of suitable decay branching ratios as well as masses and polarized production cross sections in e+e^ 
annihilation. Finally, we calculate the EDM of the r-lepton induced by the r-slepton-neutralino and 
T-sneutrino-chargino loops with complex A^, fi and Mi. 

In Section 2 we shortly review the mixing of 3rd generation sleptons in the presence of complex 
parameters. In Section 3 we give the formulae for the fermionic and bosonic decay widths of f,; and I'r- In 
Section 4 we present numerical results for the phase dependences of their branching ratios. In Section 5 
we give an estimate of the errors to be expected for the fundamental parameters and the phases of Ar , /i 
and Ml. In Section 6 we present our results for the EDM of the t. Section 7 contains a short summary. 



2 tl-tr Mixing 

We first give a short account of — ffl mixing in the case the parameters ^ and A,- are complex. The 
masses and couplings of the r-sleptons follow from the herniitian 2x2 mass matrix which in the basis 
{fL,TR) reads glEIl 



with 



1 



= A/| + (-- + sin-'evK)cos2/3 TO^ +m^, (2) 
M?^^ = M|-sin2eH/Cos2/3 m| + m2, (3) 

M|^^ = (Af?^J*==m,(A,-Ai* tan/3), (4) 

yjf = arg[^r - yu* tan /3] , (5) 

where is the mass of the r-lepton, is the weak mixing angle, tan/3 = vijvx with v\(v2) being 
the vacuum expectation value of the Higgs field Pl\{}i^\ and M^, M^,^r are the soft SUSY-breaking 

parameters of the f, system. The f mass eigenstates are (Ti,-r2) = {tl-,t'r)T^'^ with 

s**^^ cos Qr sin 6'f 
TV=\ 1 , (6) 

— sin Q=. e ^"^^ cos 0; 
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and 



cos9f = 



J\M? |2 + (to? -M? )2 |2 + (m? -M? )2 

The mass eigenvalues are 



(7) 



ml =- (A/? +Af? )T\/(A^? -a/? )2+4|A/? I2 



(8) 



The z/t- appears only in the left-state. Its mass is given by 

1 



Ml + -mz cos 2/3. 



(9) 



Eqs. (UI) and © show that the phase dependence of the mixing angle Of and the eigenvalues \ 



stems 



from the term m2 tan /? cos(iy9p + ipA^ 

and at the same time \M?^^ — M?^^\ < 
phases ipf^ and ipA^ because nir is small. 



M? I 



The phase dependence of 9f is strongest if \ Ar \ « tan/3 



The masses 



are in many cases insensitive to the 



3 Production and Decay Formulae of fi and 



The reaction e^e — s- f,;Tj proceeds via 7 and Z exchange in the s-channel. The ZfiTj couplings are 



1 



2cos9vF 
1 

2cos9vy 
1 



(cos^ Of — 2 SVC? Qw), 
(sin2 Of - 2sin2eiv), 



-e cosfc sinfc 



The reaction e'*'e 



2 cos Qw 

Vt-v-t proceeds via s-channcl Z exchange with the coupling 



2 cose 



w 



(10) 



(11) 



The cross section of e+e^ VtVt at tree level does not depend on the phases ^p^l and ^a^- The tree- 
level cross sections of the reactions €^e~ — > fif.j do not explicitly depend on the phases and <^a^, 
because the couplings C'{t* Zfi),i = 1,2, are real and in e+e^ —>■ fif2 only Z exchange contributes. The 
cross sections depend only on the mass eigenvalues ^ and on the mixing angle 9f. Therefore, they 
depend only implicitly on the phases via the cos((/9^ + fA^) dependence of mf^ ^ and 9f, Eqs. Q and (jS)). 
This holds even if one or both beams are polarized (the formulae of the cross sections including beam 
polarizations are given, e. g., in \'22\]. Of course, properly polarized e~ and e'^ beams arc a very useful 
tool to enhance some signals and reduce the background and, therefore, measure some of the observables 
with better precision [3 I28| . Information about the phases ipf^ and ipA^ separately can be obtained 
by studying the branching ratios of the fi and i>r decays into neutralinos, charginos and Higgs bosons, 
because some of them depend explicitly on the phases. It is expected that Yukawa-type corrections at 
one-loop order to the fi and i^T pair production cross sections and decay widths will not change the 
overall picture obtained in tree approximation, because they have been shown to be of the order of a few 
percent only 
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3.1 Fermionic Decay Widths of fj and v-j. 

The widths for the decays x5'''(-^t), « = 1, 2, j = 1, . . . , 4, where x° is the neutraUno and ^ ±i is 
the hehcity of the outgoing r, read 

2 / 2 2 2 \ 

r(f. - X°r(A.)) = HCT 1-^ A (12) 



with 



167rml 



+ 2(-l)^^+^i^,i^,(l4|2-|5j^r)j (13) 

where g is the weak SU{2) gauge couphng constant, k{x, y, z) — {x^ + y'^ + z'^ ~ 2xy — 2xz — 2j/z)^/^ and 
Hg = ('^li ^ ('^x" ™t)^)^, = ('^fi ~ ('^x° ~ TOt-)^) 2 . The couphngs are 



where 




with 

1 

= V2tSinewN*^. (16) 

IV = 7TiT-/(A/2mvi/ cos/3) is the r Yukawa couphng. The mixing matrices U and N are defined by 
Eqs. 1)41(1 and 151(1 in Appendices IXI and IbI Since m,- <C Wfi, we have iJ^ « iJp and, hence, to a good 
approximation, r(f, x°^(Ar)) oc (la^f) for = +^(-5) ISS- 

The width for the decay into the chargino, fi — > xj^rihj = 1,2), is obtained by the replacements 
ajj ^ij^^ij ~^ 0, m-j^o — > m^-.mr —f and A^ ^ — i in Eqs. ((T^ and ifT^ . with the couphngs ^[^ 
also given in Eqs. I(14() and ((15(1 . The width for the r-sneutrino decay Vr X^^'t is obtained by the 
replacements ajj aj", bjj ~^ 0, m^. — > m^;^, irir and Ar ^ — ^ in Eqs. p2() and ((13(1 . and that for 
the decay Dr Xj^T{XT-) by the replacements aJj ^jj^ij ^ ™n ^ "^^^ ^^i^ "^x" ^ ™x*' '^"'^^ 
couplings are now 

4 = -i=(iv,itanew-^j2), fc; = r,c/;2, ^^ = -"^^1, (17) 

with the mixing matrix V given by Eq. ((42(1 in Appendix 1X1 

As can be seen, the widths for the decays of fi and T2 into charginos and neutralinos depend on 
cos((^^ + ipA^) through mf. and df, and also on (pf, Eq. They depend also on ip^ (ip^ and (pu(i)) via 
the chargino (neutralino) masses m^-im^o) and mixing matrix U(N), see Eqs. ((40I49II (Eqs. ((501511) '). 

Xj Xj 

The widths for the I'r decays into fermions depend on the phases of (pfj, and (pu(i)- 
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3.2 Bosonic Decay Widths of T2 and I'r 

The widths for the decays of f2 and i>r into gauge bosons and Higgs bosons are given by: 

m - = \aZ^ r:2 ^ \CiKW+f2)\' (18) 

^("^ ^ - 16.m3 ^1 (^9) 



167rm? 



r(f2 ^ = ^ ^ \ciD;H+f2)\' (20) 



-r2 



r(f2 ^ il.Tl) = "3^ \C{T*Hd2)\^ (21) 

lOTTTO^^ 



r(z>. -> T^+n) = — "-'3 |C(f*W^-i>,)|2 (22) 



g^At3(TOg^,TO|^,TO^±) 

167rm?^TO^± 

9 / ? 9 9 \ 

q Kim% ,m% , m7r+ 

r(z>. ^ H+n) = ^ J3^' \c{f:H-Dr)\' (23) 



The couphngs relevant for f2 decays into the Z boson are given in Eq. (|10|1 and the couplings to the 
VF+ boson are 

C{i>;W+ fi,2) = ^(-e""^^ cos 6*^, sin 6*^). (24) 
v2 

The couplings to the Higgs bosons are more conveniently written in the weak basis {tl,te)- The 
couplings to the charged Higgs boson iJ+ are given by 

C{v*H+tl^r) = ^ (TO^tan/3-m^sin2/3,m^(tan/3|A^|e-*^'^- + |^|e'^")) (25) 

\/2mw 

The couplings C{i'*H^ti^2) of the mass eigenstates are then obtained by multiplying the couplings 
above with Tv'^ from the right. 

The couplings to the neutral Higgs bosons Hi,i = 1, 2, 3, are 

C{flH,TL) ^^(^n - i-l + sin' Qw] (cospOu - sin/JOa.) (26) 

TOvFCosp cosBw V 2 y 



C{f*„H,fR) = + ^ sin^ Qwicos f30u - sin/?02»), (27) 

niw cos p cos Idw 



C{flH,TR) = ^-—{i (sin/3|A,|e-*^^- + cos^W^-) O^, 

zmi^/cosp 

+ (|/i|e"^-02, " l^^le-'^^-^-On)}, (28) 
C(f^i?.fL) = {C{tIU,tr)\\ (29) 
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The couplings of the mass cigenstates fi are obtained by 



c•(f,*i^,f,) = 7^^• \n'\ (so) 

Oij is the real orthogonal mixing matrix in the neutral Higgs sector in the basis (0i ,(j)2,o,) ~ {\/2{TZeHi — 
vi) , \/2{JLeH2 — V2), ■\/2(sin fUmHi + cos pimH 2)) ^ where and H2 are the neutral members of the 
two Higgs doublets with hypercharge —1 and +1, respectively. Oij diagonalises the 3x3 Higgs mass 
matrix: (f)i — OijHj,i = 1,2, a = O^jHj, O'^Ad^O ~ diag{m^_^,rrijj^,rrijj^), with rriHi < 1TIH2 ^ it^h^ 
|15| . The neutral Higgs mass cigenstates Hi, i = 1,2, 3, are mixtures of the CP-even and CP-odd states, 
because of the explicit CP violation in the Higgs sector. The phase parameter ^ also introduced in 
|15[I16[[T7| does not play a role in our analysis. Therefore we put ^ = 0. 

The widths for f2 decays into the neutral Higgs bosons depend on (^^ , tpA^ and ipf and in addition on 
the mixing matrix Oij. At one-loop level Oij depends on the phases f^, ipAt and fAbi with the latter 
two being the phases of the stop and the sbottom trilinear couplings At and A^, respectively. 



4 Numerical Results 

In the following we present our numerical results showing how the fi, f2 and i>T decay branching ratios 
depend on the complex phases. In order to study the full phase dependences of the observables, we do 
not take into account the restrictions on (^^ and fu(i) from the electron and neutron EDMs. We fix the 
fi, f2 and Ur masses such that these particles can be pair produced at an e^e^ linear collider with a 
CMS energy in the range y/s = 0.5 — 1.2 TeV. Furthermore, we impose the following conditions: 

(i) m-± > 103 GcV, ttih, > 110 GeV, nif, > rn^o > 50 GeV, nif, > 80 GeV, and 

(it) \At\'^ < 3(M| -I- M|, + (m^+ -f- m| sin^ 6w) sin^ (3 — ^m^) (the approximate necessary condition 
for tree-level vacuum stability [26j ) . 

In principle, the experimental data for the rare decay 6 — > 57 lead to strong constraints on the SUSY 
and Higgs parameters in the MSSM and, in particular, in the minimal Supergravity Model (mSUGRA). 
We do not impose this constraint, because it strongly depends on the detailed properties of the squarks, 
in particular on the mixing between the squark families, which we do not take into account. 

The following parameters are necessary to specify the masses and couplings of the SUSY particles 
fi, i>r, xf and x°: M^, M^, \Ar\, (pA^, \^J■\, ififi, tan/3, M2, \Mi\, ipu{i)- Equivalently we use the 
mass eigenvalues m^^, 771^2 or the masses m^^, nio^ as input parameters instead of Mj^, M^. For the 
complete determination of the renormalization group (RG) improved MSSM Higgs sector at one-loop 
level in addition the charged Higgs boson mass to^± , the mass parameters and the trilinear couplings 
of the scalar top and scalar bottom systems Mq, Mfy, Mj^, \At\, fAt, \Ab\, fA,, and the gluino mass 
\mg\ as well as its phase Lpg = arg(mg) have to be specified ^Hj. Mixing of the CP-even and CP-odd 
neutral Higgs bosons at one-loop level is induced if A^^t and/or ^ are complex. We take to,- — 1.78 GeV, 
mt = 175 GeV, mf, = 5 GeV, mz = 91.2 GeV, sin^ Qw = 0.23, mw = mz cosQw, ct{mz) = 1/129, and 
Ois{fnz) = 0.12, where mt,b are pole masses of t and b quarks. 
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4.1 Ti Decays 

In this subsection we study the dependence of the branching ratios of fi decays into charginos and 
ncutrahnos on the phases (pA^ , <^/i and ^u{i) ■ We take mf ^ = 240 GcV. In order not to vary too many 
parameters we fix \Ar\ = 1000 GeV in Figs. 1 to 7. We assume the GUT relation |A/i| = (5/3) tan^ QwM2, 
although we take Mi complex. We focus on the decays fi — > x\ and fi Xi^t- 

We first study the i^a^ dependence of the fi decay branching ratios, because LpA^ appears only in 
the fi sector and it is the phase dependence that we are particularly interested in. In Fig. ^ we plot 
the branching ratio B{fi x\t) as a function of LpA^ for the three values mo^ = 233 GeV, 238 GeV 
and 243 GeV (corresponding to ^ 240 GeV, 245 GeV and 250 GeV), taking v?^ = (pu{i) = 0, 
1^1 = 300 GeV, tan/3 = 3, and M2 = 200 GeV. Note that B{fx Xi''') is invariant under ipA^ ~* —'^A^ 
for ip^ = {0,±7r} and (pu(i) = {0,±7r}. As can be seen, the ipA^ dependence of B{fi Xi''') is quite 
pronounced. To a large extend it is caused by a relatively strong variation of the mixing angle 0f with 
varying (pA-r- More specifically, when varying pA^ from to tt, then cos0f varies from —0.1 to —0.9 for 
mc,^ = 233 GeV, from -0.06 to -0.6 for mr,^ = 238 GeV and from -0.05 to -0.45 for m^^ = 243 GeV. 
This means that for m^^ = 238 GeV and 243 GeV fi is mainly fij-like, whereas for mp^ = 233 GeV fi 
is fi-like (ffl-like) for ipA^ ^ ^/3( ^ 7r/3). Such a strong variation of the mixing angle 6^ with pA-r can 
only occur if A/^ sa and \ Ar \ ~ tan/3, otherwise this variation is weaker. 

In the following Figs. 2 to 5 we fix nif^ — 500 GeV instead of mp^. We consider separately the two 
cases Ml < and > M^. and determine the values of and correspondingly. In Fig. [3 we 
show the tan/3 dependence of i?(fi ^ Xi''') for p^ = (solid fine), Pfj, = n/2 (dashed fine), p^ = tt (dotted 
line), with pA^ = Pu{i) = 0, M2 = 200 GeV, = 150 GeV, assuming M^ < M^. For pA, = {0,±7r} the 
branching ratios are invariant under the simultaneous sign flip {p^, fu{i)) ^ (~'/'a" ~'fiu{i))- As can be 
seen, B(fi —^ x^t) becomes almost independent of p^j, for tan/3>15. A similar behaviour is obtained for 
B{fi —> X2'^) and B(fi — > Xi'^r)- In the case of the decay fi —f Xi^r this behaviour can be understood 
by observing that the p/^^ dependence of the mass eigenvalues and the mixing matrices Uij and Vij 

changes if the value of tan/3 is changed. For the width r(fi Xi^r) we obtain from Eqs. ©, 

(dj), (HSJ and gH) 

^[1 = -e'^^{e-''^^ cos6'f cos6li - e'"^! i; sin 6*^ sin 6*1 ) . (31) 

By inspecting Eqs. (0) and (|45|l one can verify that in the limit tan/3 ^ 00 wc obtain e"*'''^ — e^'^'f' 
and e"^i 6"^^", which means that in this limit becomes independent of p^. Here note that in this 
limit Of and 9i become independent of (p^ as can be seen from Eqs. ((JJ, (jHl and (|43|1 . In the case 
of the decay into a neutralino we can see the influence of the phases p^^ and pu{i) from the approximate 
formulae 

1/1^1/^1 ™|sin2ewsin2/3cos(^^ + (p[/(i))\ 

(32) 



and 



SIT? Qw cos^ Qw 



Ml M2 



+ sin 2/3 



sin^ 61^ cos(<^^ + pu{i) ) cos^ 9;^ cos p^ 
|Mi| ^ M2 



(33) 

which hold for \M2 ± |mII ^ '^2 foi' the mass of a gaugino-like or a higgsino-like Xi, respectively. Similar 
approximation formulae hold for m^o and the mixing matrix Nij . From these formulae one can see that 
Pf_i and Pu(i) appear only in terms multiplied by sin 2/3. Therefore, in the approximation where Eqs. (|32|l 
and l|33|) hold, m^o ^ and Nij become independent of p^ and Pu(i) for large tan (3. Concerning the p^ 
dependence in general, it can be shown that m^o and Nij become independent of (^^ for tan/3 — > 00, 
because the characteristic equation of the neutralino mass eigenvalues becomes independent of p^^ in this 
limit. 
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In Figs. Ola. b wc plot the branching ratio B{fi Xit) against A/2 in the range 200 GeV < M2 < 
500 GeV for tp^ = tt (sohd hne), ipfj, = 7r/2 (dashed hne), ip^ — (dotted hne) and = — 7r/2 (dash- 
dotted Hne), taking = 0, (pu{i) = t'/'^, |mI = 150 GeV and tan/3 = 3. In Fig. |31a we assume 
Af£ < Mg, so that fi ~ fi(cos0f « — 1). This means that the couphngs are approximately |ay | ~ 
\blj\ ~ l^Ljl E^nd the decay width is essentially determined by r(fi — > x°t) oc I/^^P + In Fig.|31b 

we consider the case > Mj^. In this case we have \a\j\ — l^^yl — and ~ Yr\Uj2\- 

This means that the decay fi xT^t suppressed, because now fi ~ ffl; {cos0f w 0) and the tiXi^t 
coupling is nearly proportional to the small Yukawa coupling 1^. Therefore, B{fi Xit) in Fig. 
is larger than in Fig. |Sla. In both cases there is a significant variation with ipf^ . The (p^ dependence 
of B{fi — > Xi''') in Figs.|31a, b is caused by an interplay between the dependence of the mass and 
mixing character of the fi and that of the Xi- The M2 dependence can be understood by noting that 
for M2 ~ 200 GeV the lightest neutralino has a sizable gaugino content, which decreases for increasing 
Af2- For our parameter choice Xi becomes mainly higgsino-like for Af2^300 GeV. Near A/2 ~ 440 GeV 
the decays into gaugino-like neutralinos become kinematically forbidden, which causes the increase of 
B{fi Xit) for A/2>400 GeV. 

We have studied the 1^9^ dependence of B{fi Xi"^) a^lso for other values of and have found that 
it is less pronounced if |/i|>A/2 and that it is stronger if « M2 or |/i|<|A/i|. As shown in Fig. |21it is 
stronger for low tan /?. 

In Figs. 01a, b we show the (pu{i) dependence of B{fi Xi''') foi' l/^l ~ 150 GeV, tan/3 = 3 and 
(pA^ = 0, for Pf^ = TT (solid line), p^ ~ tt/2 (dashed line), (/9^ = (dotted line) and p^ = —tt/2 
(dash-dotted hne). In Fig. jUa we take A/^ < and A/2 = 280 GeV. Fig. Hb is for A/^ > A/^ and 
A/2 = 380 GeV. Although the pu(^i) dependence of B{fi Xi''') stems only from the Pu{i) dependence 
of the Xi 2 parameters, it is quite pronounced. It is essentially explained by the dependences of A'^n 

and A^i2, which enter in the couplings f£-^ and fj^^ (see Eqs. (|14|) - For example, the minimum 

of B{fi Xit) in Fig. ^b at pu{i) ~ 37r/4 (— 37r/4) for p^^ = 7r/2 {—tt/2) is caused by a corresponding 
minimum of |iVii|. 

We have also studied how the branching ratios B{fi X2 3T) and B{fi Xi^t) vary as functions of 
the phases. As an example wc show in Fig.jSlthese branching ratios as functions of p^ for Pij(i) = pa^ = 0, 
A/2 = 280 GeV, = 150 GeV and tan/3 = 3, assuming A/^ < A/^. For this set of parameters all 
branching ratios shown have a significant p^^ dependence. Their behaviour can be understood in the 
following way: If we first consider B{fi Xi^t) oc the p^ dependence of \Uii\ follows from 



where Oi is the mixing angle of the chargino mixing matrix Uij defined in Eq. (|43|l . The mass squared 
difference m?+ — m-+ decreases for p^^ — > tt, which can be seen from Eq. I^?^. therefore, also |J7ii| 

decreases. The behaviour of i3(fi Xi,2,3''') can be understood by noting that Xi,2,3 have large higgsino- 
components. Varying p)^ from to tt essentially interchanges the Hi and //^ components of x? 2 3- This 
causes the variation in the branching ratios, because fi couples to the component of Xi* but not to 
the H2 component. 

It is expected that p^ and p>u{i) will be determined by measuring suitable observables of the chargino 
and neutralino sectors |27]. The and Pu{i) dependences of the various fi decay branching ratios, 
however, will give useful additional information for the precise determination of (/s^ and Pu(i) thereby 
provide further tests of the MSSM with complex parameters. This may also be helpful for resolving the 
ambiguities encountered in the studies about the parameter determination of the chargino and neutralino 
sectors |77| . 

An additional observable which is very sensitive to the SUSY parameters of the n and Xfe systems is 




(34) 
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the longitudinal polarization of the outgoing r-lcpton in the decays fi x'j''' For the fi decays into 
neutralinos it is defined as 

where the last equation holds in the limit m,- ^0. R,L denote At = +^1^57 respectively. 

We show in Figs. IHla, b the longitudinal polarization of the r in the decays fi x^t and fi 
X2''', respectively, as a function of (pA^ for m^;^ = 233 GeV (solid line), 238 GeV (dashed line) and 
243 GeV (dotted hue), which correspond to Mi = 240 GeV, 245 GeV and 250 GeV, respectively. The 
other parameters are M2 = 200 GeV, — 300 GeV, tan/3 — 3, ipfj, — (pu{i) = 0. The behaviour of 
VriXi''') in Fig. Ela follows from the change of the mixing angle 9f with varying (pA^, as described in 
the discussion of Fig. The behaviour of VrixHT) in Fig. IHIb can be understood by noting that in this 
case X2 is mainly a which couples only to the tl component of fi and that this component strongly 
increases for LpA^ ^ tt as can be seen from Eqs. © and (U)). 

In Fig. [Tlwe show the longitudinal r polarization in the decays fi x^t and fi x^'t ^ function 
of (fifj^. Here we have taken rrif^ — 500 GeV and the other parameters M2 = 350 GeV, \fj,\ = 150 GeV, 
tan/? = 3, ipu{i) — — 0. As we have chosen AI^ < AIj^, ti is mainly a and Vr is negative for 
ip^ > 37r/10 due to the very small r Yukawa coupling. For (^^ 0, the fLTLXi.2 couplings \f£i\ and 
1 7^2 1 decrease monotonically, because Xi 2 mainly higgsino-like and changing the phase (^^ from tt 
to implies essentially a decrease of their gaugino components as well as exchanging the component 
with the H2 component. This leads to a change of the sign of Vr- I/L2I has a maximum at (/j^; w 37r/4, 
which is clearly seen in the minimum of Vrix^T) sa —0.6 for this value of ip^. 



4.2 T2 Decays 

As we have seen in the previous subsection, the branching ratios for the fermionic fi decays depend on 
the phase ipA^ only via the cos((/?^^ + ip^) dependence of the mass m^^ and the mixing angle Of. We 
consider now the bosonic T2 decays where the couplings to the Higgs bosons explicitely depend on the 
phases pA^ and (see Eqs. to (I^Ull V The decay widths into W^, Z and Higgs bosons are enhanced 
by choosing and/or \Ar\ large |28| . 

As already mentioned, the RG improved Higgs sector is determined by the parameters mfj± , tan /?, 
\At\, \Abl pf,, PA,, ipA„ Mq, Mjj, Mf,, \mg\, pg, |Mi|, puil) and M2 M We fix Mq = Mfj = 
Mfj = MsusY- The amount of the CP violating scalar-pseudoscalar transition in the neutral Higgs mass 
matrix is proportional to the parameter 

'^^^ = 1287r2m^M|^,^^^"(^''+'^^^)' ^^^^ 

where f = t,b |15l llHj. This means that significant CP violating effects in the Higgs sector can be 
expected if \Af \ > Msusy and | sin((y5^t + p>Aj )| ~ 1- As we focus on the pA^ and the p^ dependence 
of the observables, we fix the phases p)At = pjg = pAt = and we take \At\ ~ \Ai,\ — 800 GeV, 
Msusy = 600 GeV, \m~g\ = ias{\mg\)/a2)M2 (with a,{Q) = 12V((33 - 2n/) ln(QVA', )), «/ being the 
number of quark fiavors). For this choice of parameters mixing between the CP-even and CP-odd Higgs 
bosons at one loop level occurs only if p^ ^ {0,±7r}. Therefore, we can control the influence of explicit 
CP violation in the Higgs sector with the parameter . With this choice of parameters the constraint 
from the p-parametcr on the t and b masses and mixings, Sp{t — b) < 0.0012, is always fulfilled |29j . 

For large tan/3 the allowed range of |^| is restricted by the two-loop contributions to the EDMs 
of electron and neutron PO]. For example, for tan/3 = 40, = n/2, mjj±<200 GeV and the other 
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parameters as fixed above the EDMs give the restriction |/i|<600 GeV. Therefore, we also fix |^| — 
600 GeV. 

Ill tlie following we give some numerical examples which show the dependence of the branching ratios 
for f2 ^ TiHi,i = 1,2,3, on i^a^, tan/3 and mjj±. We take \Ar\ = 900 GeV, Ma = 450 GeV and 
fuii) = 0- We consider the case > M^, where fa is mainly fi,-like and fi is mainly fu-^ike. In this 
case the decays fa — > W~Vr and fa H~Vt are kinematically forbidden. 

In Figs.lSla, b we show the branching ratios for various fermionic and bosonic fa decays as a function 
of ipA^ for (/J^ = and 7r/2, taking tan/3 = 30, TO_f/± = 160 GeV, rrif^ = 240 GeV, rrif^ = 500 GeV and 
the other parameters as specified above. As can be seen, the branching ratios of the decays fa — > ^fi,2,3'?i 
show a pronounced {^pA^^^^J) dependence. The behaviour of these branching ratios can be understood 
by examining the approximate formula for the coupling squared for fa fi Hi, 

|C(fa*i/,fi)P ~ |C(f£i/,ffl)p (1 - 2 sin2 9^ cos^ 0^(1 + cos 2(arg[C(f£i^,fi^)] + iff))) (37) 

with 

|C(f2iI,fH)P 0, lY,'(^{\^,\'-\Ar\')0l + \Ar\'-2MA^\02. 

X (Ou cos{(pA^ + (p^) - Osi sm{LpA, + (p^)j ^ , (38) 

which follows from Eq. (|28|l and H30() . Here we have omitted terms proportional to (C(f£i/ifL) — 
Cij'^Hifii)) and cos^. Eqs. and show that a significant phase dependence of the fa fi-ffj 
branching ratios can be expected for large tan (3. Moreover, also the dependence of the Higgs mixing 
matrix elements influences in a significant way the behaviour of i3(f2 — > fi Hi). For (p^ ~ 0, for 
example, we obtain On w -0.262, Oai ~ -0.965, ©31 = Cia = O22 = 0, ©32 1, C13 ~ 0.965, 
C23 w -0.262, O33 = 0, mn, = 115.74 GeV, m^, = 138.48 GeV, m^, = 139.14 GeV. The ipA^ de- 
pendence of i3(fa — > HiTi) follows essentially from the cos((^^^ + f^) term and the first two terms of 
Eq. (|38|) . The minimum of i3(f2 — > -ffifi) at pA^ = (Fig. |Hla) follows from a partial cancellation of 
the terms in Eq. (|38|l (or, equivalently, from a partial cancellation of the last two terms of Eq. H28|) . 
see also Fig. below). The cos{ipA^ + Ptj) term and the first two terms of Eq. H38|l determine also the 
ifA^r behaviour of i3(fa H^ti). The pA^ dependence of i3(fa H2T1) follows from the last factor 
of Eq. ^ and the first term of Eq. As for Fig. |Hlb, for = 7r/2 we obtain On « -0.106, 

O21 ~ -0.992, w 0.066, da ~ -0.230, O22 ~ -0.040, 032 ~ -0.972, dg « 0.967, ©as « -0.118, 
C33 ~ -0.224, m.H, = 117.09 GeV, tuh^ = 138.48 GeV, = 139.14 GeV. The pA^ dependence of 

i3(fa — * Hifi) is now different from that in Fig. |Sla. In the case of i3(fa i?ifi) the cos(<pAt + fti) 
term becomes —suvipAr a-nd it is multiplied by a much smaller factor, which explains the relatively flat 
LpA^ dependence. The behaviour of i3(fa — > i?2'?i) and i3(f2 ^ H^fi) can be explained in an analogous 
way. For comparison we also plotted the branching ratios of fa Zf\ and of some of the decays into 
charginos and neutralinos. The ^pA^ dependence of fa —* Zt\ essentially drops out (see Eq. l(Tn|) ') and 
that of the fermionic decays disappears due to the large value of tan /3 for which is insensitive to (pA^ ■ 

We also studied the tan /3 dependence and the m^± dependence of the fa decay branching ratios into 
neutral Higgs particles. For tan/3 these branching ratios vanish (cx tan/3), whereas for tan/3 > 10 
they depend only weakly on tan/3. The m^± dependence of the branching ratio -B(fa H\f\) is shown 
in Fig.l^lfor LpA^ — 0, 7r/2, tt. At m^± = 150 GeV and pA^ ~ this branching ratio practically vanishes. 
The reason is that the coupling C(f£-ffif/f) practically vanishes for this set of parameters due to a 
cancellation of the last two terms in Eq. At this point also a level crossing of H\ and H2 occurs. 

We see that this branching ratio is sensitive to p}A^ for m/^± < 250 GeV. 
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4.3 Ur Decays 

The decay widths for Vt decays into charginos and neutrahnos are independent of ipA^ ■ The decay widths 
for v-r — > XfeT depend on ip^, those for Vt — > Xfc^r depend also on '■Pu(i)- We first assume < Af^. 
which leads to a sneutrino mass m^^ ~ 229 GeV for irif^ — 240 GeV, rrif^ = 500 GeV and tan/3 — 3. In 
this case the decays i^r — >■ W'^fi and — > H^fi are kinematically forbidden. 

We show in Figs.llOla and b the branching ratios for the decays into Xi^^t, X2^t and x^r as functions 
of ifi^ and ipu{i), respectively, for Ms = 500 GeV, = 150 GeV, tan/3 = 3, and \Ar\ = 1000 GeV. 
In Fig. II Ol a we take ifu{i) — and in Fig. llOl b wc take ip^^ = 0. As can be seen, the branching ratio 
for Dt Xi^T decreases for tt, whereas those for Dj- -~> X2^t and Vr Xi''' increase. The 

decay widths T{vt Xi^^t) and T{vr X2^t) decrease for (^^ — > tt, because the matrix elements |iVi2| 
and |iV22| decrease for ipf^ tt. The matrix element |Vii| entering the decay width r{i>T xt''') ^1^° 
decreases, see Eqs. (|42|l and (|44|) . However, as r{i/r Xi^''") and T{i>T X2^t) decrease more slowly 
than the total decay width, the corresponding branching ratios increase for ip^ —y tt. In Fig. llOl b the 
branching ratio B{Dr Xi^-r) decreases for (/^[/(i) tt and i?(i>T ~^ Xi't) increases. The reason is that 
|7Vii t&nQw ~ Afi2| and hence the width r{i'r — > Xi'^t) rapidly decreases for (pu{i) ^ -B(i>t- Xi't) 
increases due to the decrease of the total decay width. 

In the case M^ — A'I i^>mw, 'mu+ also the bosonic decays Vr — ^ W~^fi, if+fi are kinematically allowed. 
Consequently, the branching ratios of the fermionic decays are reduced. It turns out that in most 
cases the bosonic decay widths are almost independent of the phases; only in the region ^ mci^ < 
|Mi_2| a significant dependence on the phases is possible. For small tan/3 the phase dependence of the 
width r{i/r H~^fi) tends to be suppressed, because of the small Yukawa coupling, sec Eq. H25f) . For 
large tan/3 the term nir A* ta.n [3 in Eq. (|25|l dominates and T{i>r H^fi) cx \sm6fC{i>*H'^Tft)\'^ oc 
I sin0fmT-yl* tan/3p becomes essentially independent of the phases. Note here that Of is hardly sensitive 
to the phases because M? — M|, ^ mr\Ar — fi* tan/3| in these scenarios. The phase dependence of the 
width r{Dr W^Ti) is caused only by the phase dependence of cos6'f (see Eq. (|23J) and is again weak 
by the same reasoning as above. 



5 Parameter Determination 

We now study the extent to which one can extract the underlying parameters from measured masses, 
branching ratios and cross sections. In the following we assume that an integrated luminosity of 2 ab~^ 
is available. At a high luminosity collider like TESLA one can expect that this amount of integrated 
luminosity will be accumulated in four years of running Our strategy is as follows: 

1. Take a specific set of values of the MSSM parameters. 

2. Calculate the masses of fi, xf, the production cross sections for e+e^ fifj and branching 
ratios of the fi decays. 

3. Regard these calculated values as real experimental data with definite errors. 

4. Determine the underlying MSSM parameters and their errors from the "experimental data" by a 
fit. 

We have checked that inclusion of the data on the mass, production and decays of i>r does not further 
improve the accuracy of the underlying parameters to be determined. The reason is that the expected 
relative errors of the data in the sneutrino sector are larger than those in the stau sector |31l • 
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Table 1: Calculated masses and their assumed errors (in GeV). 



tan /? = 3 


tan/3 = 30 


TUf, = 155.0 ±0.7 
m^o = 125.6 ±0.17 
m^o = 253.5 ± 0.24 
m^l = 194.0 ±0.06 


TUf^ = 352.6 ± 1.2 
m^o = 205.6 ±0.11 
m^o = 343.1 ±0.51 
m~l = 340.9 ±0.25 

X5 


rrif, = 150.6 ±2.1 
m^o = 133.2 ±0.56 
m^o = 258.0 ±0.73 
m^l = 210.0 ±0.19 


= 355.7 ± 3.6 
m^o = 214.3 ±0.35 
TO^'o = 331.4 ± 1.4 
m-+*'= 331.6 ±0.72 

X9 



Table 2: Branching ratios of f2 decays calculated for = 150 GeV, = 350 GeV, A^- = -800 i GeV, 
M2 = 280 GeV, fi = 250 GeV and ipu{i) = 0. We show only branching ratios larger than 10^"^. 



tan/3 






















3 


0.116 


0.423 


0.001 


0.002 


0.438 


0.008 


0.002 


0.008 


0.003 





30 


0.107 


0.195 


0.036 


0.008 


0.135 


0.019 


0.044 


0.393 


0.062 


0.001 



We have taken the following input parameters for the calculation of these observables: Mj^ = 150 GeV, 
Ml = 350 GeV, Ar = -800 i GeV, M2 = 280 GeV, /x = 250 GeV and ipu(i) = 0. We have considered the 
cases tan/3 = 3 and 30. The Higgs sector has been fixed with mH+ — 170 GeV (160), m^o = 151.4 GeV 
(138.5), m,,o = 113.3 GeV (115.7), toj^o = 155.6 GeV (139.1) and sin a = 0.432 (-0.26) in case of 
tan/3 = 3 (30). Here /i", A^, a are the lighter CP-even Higgs boson, the heavier CP-even Higgs 
boson, the CP-odd Higgs boson and the mixing angle of the CP-even Higgs bosons, respectively. Here 
we focus on the determination of the phase tpA^ of A^-, therefore, we neglect mixing of the CP-even and 
CP-odd Higgs states. We have taken the relative errors of stau masses, chargino and neutralino masses 
from 01221 J which we rescale according to our scenario; in case of tan /3 = 30 we have taken into account 
an additional factor of 3 for the errors (relatively to tan /3 = 3) due to the reduced efficiency in case of 
multi T final states as indicated by the studies in [23 . We take the errors of the Higgs mass parameters 
as Am/jO = 50 MeV, Amjyo = AniAo = Am//+ = 1.5 GeV (HI for tan/3 = 3 and 30. For the branching 
ratios and the production cross sections we have taken the statistical errors only. We give the values of 
the calculated masses and assumed errors in Tableland those of the calculated branching ratios of f2 
decays in Tablc|5J fi decays only into rx? for both values of tan/3, because this is the only channel open. 

For the determination of the stau parameters we have used the information obtained from the mea- 
surement of the stau masses at threshold and the production cross sections of fifj pairs at y/s = 800 GeV 
for two different (e~,e+) beam polarizations (P_,P+) = (0.8,-0.6) and {P-,P+) = (—0.8,0.6). Here 
we have assumed that a total effective luminosity of 250 fb""'^ is avaible for each choice of polarization. 
The cross section measurements are important for the determination of |cos0fp as can be seen from 
Eq. (|10|) and the formulae for the cross sections in 0. In addition we have used the information from 
all branching ratios in Table |3 (with corresponding statistical errors) . These branching ratios together 
with the masses and cross sections form an over-constraining system of observables for the parameters 
M£, Mj^, ^eAr, QmAr, ^efj,, 3m/i, tan/3, 3?eMi, QmMi, M2. We have determined these parameters 
and their errors from the "experimental data" on these observables by a least-square fit. The results 
obtained are shown in Table |2| As one can see, all parameters can be determined rather precisely, tan/3 
can be determined with an accuracy of about 2% in the case of tan/3 = 30 and about 1% in the case 
of tan/3 = 3. The relative error of the remaining parameters except A^ is about 1%. For Ar we obtain 
the errors AQmAr/\Ar\ « 9%, A'^eAr/\Ar\ « 22% in the case tan/3 = 3, and A'^mAr/\Ar\ « 3%, 
A3?eAr/|Ai-| « 7% in the case tan/3 = 30. At first glance it might be surprising that the errors of the 
stau parameters are relatively small in case of large tan /3, despite the fact that the assumed errors of 
the masses are larger for large tan/3. The error of Ar even decreases. The reason for this is the large 
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branching ratio for f2 A'^fi in the case tan/? — 30 and the input parameters chosen (see Table EJ, 
which gives a strong constraint on \Ar\. For the determination of At it is important that the f2 decays 
into neutral Higgs bosons are kinematically allowed, because their couplings to the staus are practically 
proportional to At tan /3. Otherwise one would have to include the decays of the heavier Higgs bosons to 
get additional information on At from their decays into staus. This will be discussed in a forthcoming 
paper PS]. Additional information could also be obtained at a collider. In case of sizable CP 

violating phases fifi pairs can be produced at the resonances of both heavier neutral Higgs states H2,3 
[TO] whereas in case of CP conservation fifi pairs can only be produced at the resonance but not at 
the A^ resonance |36| . 

In the procedure described above we have determined the errors of the fundamental parameters 
assuming an integrated luminosity of 2 ab~^, taking the expected experimental errors of the masses from 
the Monte Carlo studies in [H| and rescaling them to our scenario. It is clear that further detailed 
Monte Carlo studies including experimental cuts and detector simulation are necessary to determine more 
accurately the expected experimental errors of the observables for our scenario, in particular the errors 
of the stau decay branching ratios. Such a study is, however, beyond the scope of this paper. Instead 
we have studied how our results for the errors of the fundamental parameters are changed when the 
experimental errors of the various observables are changed: we have redone the procedure doubling the 
errors of the masses and/or branching ratios and/or cross sections. Clearly we have found that the errors 
of all parameters are approximately doubled if all experimental errors are doubled. Moreover, in this 
way we can see to which observables an individual parameter is most sensitive. Concentrating on the 
stau sector we find that the precision of M|, and M| is sensitive to the stau mass determination at the 
threshold as well as to the measurement of the total cross sections in the continuum. The accuracy of At 
is most sensitive to precise measurements of the branching ratios, especially to those for the decays into 
Higgs bosons. The precision of /i is more sensitive to the errors of chargino and neutralino masses than to 
the errors of the stau observables. In the case of large tan/3, the precision of tan/3 depends significantly 
on the precision of the stau cross sections and to a lesser extent also on that of the stau decay branching 
ratios. 

In our procedure we have also determined the expected errors of 5Re/i, 3to^, tan/3, SReAfi, SjtoMi, 
M2 using also the information obtainable from mass measurements of charginos and neutralinos. As one 
can see in Table|31 the results are quite satisfactory. Once these parameters together with the Higgs mass 
and mixing parameters are precisely determined in the chargino, neutralino and Higgs sectors, one can 
then include them as input values in the determination of the parameters of the stau sector. This will 
in turn improve the accuracy in the determination of 5Re(^T) and 3771(^47-). Note that this accuracy of 
the paramters at the weak scale allows also a rather precise determination of parameters at a high scale, 
e.g. the GUT scale, and hence the reconstruction of the parameters of an underlying theory at this high 
scale |S7). 



6 Electric Dipole Moment of the r-lepton 

The MSSM with complex parameters implies also a possible electric dipole moment (EDM) of the r- 
lepton, which is induced by chargino-sneutrino as well as stau-neutralino loops. For the calculation of 
the T EDM we use the corresponding formulae given in |12j for the electron EDM by replacing me by to^ . 
It turns out that the natural range for the t EDM is O(10~^^) — O(10~^^) ecm. This is demonstrated 
in Fig. Illl where we show the r EDM (It corresponding to some of the scenarios discussed above. This is 
about 5-6 orders of magnitude below the current experimental limit: |c?^^p| < 3.1 • 10~^^ ecm |38) . 

The dominant contribution stems from the chargino loops as in case of electrons. However, for the 
r EDM the neutralino loop is much more important than in case of the electron due to the fact that 
uIt ^ TOe- Its modulus can reach about 10% of the chargino-loop contributions as can be seen in Figs.lTTb 
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Tabic 3: Extracted parameters from the "experimental data" of the masses, production cross sections 
and decay branching ratios of fi. The original parameter point is specified by: — 150 GeV, = 
350 GeV, Ar = -800 i GeV, A'h = 280 GeV, /i ^ 250 GeV and = 0. 



tan/3 


3 


30 


M| [GeV^] 
m| [GeV2] 


2.25 -lO-^ib 2.2 -10^ 


2.25 -lO^i 6.0 -10^ 


1.225 -lO^i 4.3 -lO^ 


1.229 -lO^ib 7.0 -lO^ 


^e{Ar) [GeV] 


-8.0 ± 180 


8.0 ± 55 


Qm{Ar) [GeV] 


-800 ± 70 


-800 ± 21 


3?e(/i) [GeV] 


249.9 ± 0.26 


249.9 ± 0.6 


3TO(Ai) [GeV] 


2.4 ± 1.7 


-0.2 ± 3.8 


tan (3 


2.999 ± 2.7 -10-2 


29.9 ± 0.70 


JRe(Mi) [GeV] 


140.9 ± 0.21 


140.6 ± 0.63 


3to(Mi) [GeV] 


-0.7 ± 3.4 


0.16 ± 1.0 


M2 [GeV] 


280 ± 0.29 


280 ± 1.0 



and f. The solid line shows the total r EDM, the dashed line the chargino-loop contributions and the 
dotted line the neutralino-loop contributions. In the other plots of Fig. the r EDM is identical 
to the neutralino-loop contributions, because in these scenarios (p^ ~ and hence the chargino-loop 
contribution vanishes. 



7 Summary 

In this paper we have presented a phenomenological study of r-sleptons 7^ and r-sneutrinos in the 
Minimal Supersymmetric Standard Model with complex parameters Ar, fi and Mi. We have taken 
into account explicit CP violation in the Higgs sector induced by ti and hi loops with complex and 
complex trilinear coupling parameters At and Ab- We have analysed production and decays of the fi 
and at a future e'^e" linear collider. We have presented numerical predictions for the fermionic and 
bosonic decays of fi, f2 and Vt- We have analyzed their SUSY parameter dependence, paying particular 
attention to their dependence on the phases fA^, f^i and </5£/(i)- For tan/3<10 the phase dependence of 
the branching ratios of the fermionic decays of fi and Dr is significant whereas it becomes less pronounced 
for tan/3 > 10. The branching ratios of the f2 decays into Higgs bosons depend very sensitively on the 
phases if tan/3>10. Quite generally one can say that the decay pattern of the fi and D^- becomes even 
more involved if the parameters A^-, /i and Mi are complex and if mixing of the CP-even and CP-odd 
Higgs bosons is taken into account. 

We have also given an estimate of the expected accuracy in the determination of the MSSM parameters 
of the fj sector by measurements of the masses, branching ratios and cross sections. We have considered 
the cases tan/3 = 3 and tan/3 = 30. We have found that on favorable conditions the accuracy of the 
parameter A^ can be expected to be of the order of 10% and that of the remaining stau parameters 
in the range of approximately 1% to 3%, assuming an integrated luminosity of 2 ab^^. In addition we 
have considered the electric dipole moment of the r-lepton induced by the complex parameters in the 
stau sector as well as the chargino and neutralino sectors. We find that it is well below the current 
experimental limit. 
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A Chargino Masses and Mixing 

The chargino mass matrix in the weak basis is given by 121j 



AI2 is the SU{2) gaugino mass parameter, cp and sp arc shorthand notations for cos/? and sin/3, respec- 
tively. This complex 2x2 matrix is diagonalized by the unitary 2x2 matrices U and V: 




(39) 



U*McV^ = diag{m-±,m^±), 



< m-± < m-± . 



(40) 



The unitary matrices U and V can be parameterized in the following way: 




(41) 




(42) 



with 



tan 26* 



2V2mw[Micl + |Af|2s2 +M2|Ai|sin2/3cosv3^]i/2 



(43) 



tan 26*2 



M| - |^|2 - 2m2^cos2/3 
2V2mw[Misl + l/ipc^ + Af2|M|sin2/3cos(^,Ji/2 



(44) 



M| - |^|2 + 2m2^cos2/3 




(45) 




(46) 




(47) 
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tan 72 = sin cos Lp^^ 



iVhrri^ sin 2/3 
|Ai|(m2 -M|) 



(48) 



where — 7r/2 < 6*1^2 < 0. The mass eigenvalues squared are 



Xl,2 



= - M. 



2ml, T - ImH' + 4m^ cos^ 2/? + Aml^{Ml + |//|^ 



■ 2A'/2|mI sin 2/3 cos v?^) 



(49) 



B Neutralino Masses and Mixing 



The ncutrahno mass matrix in the weak basis {B, , Hf, H2) is given as [Tl 121): 

/ iMile^'^^'d) -mzswcp mzswsp \ 

A/2 mzcwcf} -mzcwsp 

—■mzswC/3 mzcwcp —\fi\e'-'^'' 

\ mzswsp -mzcwsp -\^\e"^'^' / 



Mn 



(50) 



where Mi is U{1) gaugino mass parameter, with (pu{i) being the phase of Mi; cw and sw are shorthand 
notations for cosQw and sinGiv, respectively. This symmetric complex mass matrix is diagonalized by 
the unitary 4x4 matrix N: 



N*A4nN^ = diag{m^o, . . . ,m^o), 



< m^o < . . . < m^o 

— Ai — — A4 



(51) 
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Figure 1: Branching ratio of fi — > xjr as a function of ipA-r for rrif-^ = 240 GeV, mp^ = 
233 GeV (solid line), 238 GeV (dashed hne), 243 GeV (dotted hne), and = ifu{i) = 0, = 300 GeV, 
l^^l = 1000 GeV, tan/3 = 3, and M2 = 200 GeV. 




Figure 2: Branching ratio of fi Xi^ as a function of tan/3 for m^^ = 240 GeV, rrif^ = 500 GeV, 
(p^ = O(solid line), 7r/2(dashed line), 7r(dotted line), with the other parameters ifA^ = 'fu(i) = 0, 
A/2 = 200 GeV, \n\ = 150 GeV, and \Ar\ = 1000 GeV, assuming < Af^. 
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Figure 3: Branching ratio of fi Xi^ a-s a function of M2 for = 7r(solid line). 7r/2(dashcd line), 
O(dotted line), -7r/2(dashdottcd line), ipA^ = 0, (pu{i) = 7^/2, ruf^ = 240 GeV, m^^ = 500 GeV, 
|/x| = 150 GeV, tan/3 = 3, and \Ar\ = 1000 GeV, assuming a) Mi < Af^, b) Mi > M^. 




Figure 4: Branching ratio of fi Xi''' ^ function of ^m^i) for ip^ = 7r(solid line), 7r/2(dashed line), 
O(dotted line), -7r/2(dashdotted line), (pA^ = 0, rrif, = 240 GeV, rrif^ = 500 GeV, \^j,\ = 150 GeV, 
tan/3 = 3, and \ Ar \ = 1000 GeV, assuming a) Mi < Af^, M2 ^ 280 GeV, b) Mi > Af^, A/2 = 380 GeV. 
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Figure 5: Branching ratios of fi — * Xi 2 3''' ^'^d ^1 ~* Xi as a function of for (/';7(i) = Va^. = 0, 
TOf, = 240 GcV, = 500 GeV, M2 '= 280 GeV, = 150 GeV, tan/3 = 3, and = lOOo'^GeV, 
assuming A/^ < Afjj. 




Figure 6: Longitudinal r polarization, defined in Eq. I|35|) . for a) fi ^ Xi''' a-nd b) fi — > X2''' ^is a function 
of The parameters are ruf^ = 240 GeV, ?7ip^ = 233 GeV (solid fine), 238 GeV (dashed fine), 243 GeV 
(dotted line), = = 0, M2 = 200 GeV, \n\ = 300 GeV, tan/3 = 3, and \Ar\ = 1000 GeV. 
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Figure 7: Longitudinal r polarization, defined in Eq. H35|) . for fi — > x^''' s-nd "^i ~^ X2''" function of (p^. 
The parameters are ipu{i) = </5a, = 0, ruf, = 240 GcV, m^^ = 500 GeV, M2 = 350 GeV, = 150 GeV, 
tan/3 = 3, and = 1000 GcV, assuming < M^. 




Figure 8: Branching ratios of f2 Hi,2,3T^i, "^2 ^ Zti, T2 — > Xi 2''" ^^'^ ''^2 ^ Xi as a function of 
ifiA^ for a) iy9p = and b) = 7r/2, with the other parameters mf^ = 240 GeV, rrif.^ = 500 GcV, 
mj± = 160 GeV, l^f] = 600 GeV, Ivh = 450 GeV, ipu{i) = 0, tan/3 = 30, and = 900 GeV, assuming 
Ml > M^. 
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Figure 9: Branching ratio of f2 Hifi, as a function of mjLf± for ipA^ = O(solid line), 7r/2(dashed line), 
7r(dotted line), (pf, = (pui^i) = 0, m^, = 240 GeV, m^, = 500 GeV, Vl = 600 GeV, \Ar\ = 900 GeV, 
tan/3 = 30, and A/2 = 450 GeV, assuming > Af^. In the grey area the condition rriH^ > 110 GeV is 
not fulfilled. 





Figure 10: Branching ratios of Xi ^-nd ^ ^ ^ function of a) ip^ for <pu(\) = a-^d 

b) P^u{\) for '/'m = ^- The other parameters are ra^^ = 240 GeV, 771^3 = 500 GeV, A/2 = 500 GeV, 
= 150 GeV, tan/3 3, \Ar \ = 1000 GeV and ipA^ = 0, assuming A/^ < M^. 
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Figure 11: dr (in 10~^"'^ecm) corresponding to a) Fig. 1, b) Fig. 5, c) Fig. 4a with tp^ = 0, d) Fig. 4b 
with Lp^ ~ 0; e) Fig. 8a, and f) Fig. 8b. The hncs in a) correspond to = 233 GeV (soUd hne). 
238 GeV (dashed hne), 243 GeV (dotted hne). The hncs in b) and f) correspond to: total r EDM (sohd 
hne), chargino-loop contribution (dashed hne) and ncutrahno-loop contribution (dotted hne). 
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